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A.  Project  Objectives 

The  goal  of  this  research  project  is  to  make  a  source  of  coherent  THz  radiation  using  Difference 
Frequency  Generation  (DFG)  technique  by  Photomixing  of  two  optical  beams  with  close 
frequencies  in  a  thin  film  nonlinear  layer.  In  Phase  I  the  main  goal  is  to  formulate  a  design  of 
such  a  system  which  would  provide  efficient  energy  transfer  from  optical  laser  beams  into  THz 
wave. 

Thin  film  layer  is  a  part  of  the  multi-layer  waveguide  that  supports  THz  mode.  Design  of  the 
multilayer  waveguide  must  provide  Phase  Matching  Condition  (PMC)  between  THz  and  optical 
waves.  If  this  condition  satisfied,  optical  waves  will  efficiently  transfer  their  energy  into  THz 
wave. 

In  this  research  project  we  consider  the  nonlinear  medium,  where  two  main  waves  with 
frequencies  g>i  and  0)2  are  mixed,  to  possess  a  second-order  nonlinearity.  The  nonlinear  part  of 
material  polarization  is  expressed  as  PiNL  =  2dykEjEk,  where  dyk  is  the  second-order  optical 
susceptibility  tensor  (i,j,k=x,y,z),  and  Ej  is  the  electric  field.  In  general,  we  can  have  all  four 
possible  output  frequencies,  namely:  2t0i ,  2(02  ,  (Oi  +  (02  ,  and  toi  -  (02  (DFG).  At  the  end  of  this 
Phase  I  project  we  clarify  this  possibilities  and  also  formulate  a  way  which  will  lead  to  the 
waveguide  design  for  Double  Harmonic  (DH)  case  0Joutput=2(Oi.  This  will  open  the  way  for  much 
wanted  efficient  UV  generation  in  200-300nm  spectrum  region. 

Confining  two  idler  optical  beams  into  a  nonlinear  thin  film  waveguide  will  definitely  induce  the 
maximum  nonlinear  polarization  in  the  medium.  Having  correct  multi-layer  waveguide  design 
should  lead  to  efficient  excitation  of  the  DFG  THz  mode  by  this  nonlinear  polarization.  It  is  also 
important  to  note  that  a  planar  design  of  efficient  THz  source,  which  we  seek  in  this  project,  will 

have  an  ideal  case  for  future  integration  of  this  source  with  planar  THz  sensing  platforms 
and  general  THz  electronics. 

In  Phase  I  of  this  SBIR  project  we  have  addressed  the  following  technical  objectives: 

1.  Design  a  multi-layer  waveguide  comprising  of  the  optical  and  THz  wave  confinement 
layers 

2.  Design  a  grating  coupler  embedded  in  the  thin  film  optical  layer.  Determine  what  lens 
configuration  provides  the  optimize  focusing  of  the  idler  beams. 

3.  Predict  characteristics  of  the  future  THz  source.  Theoretically  estimate  efficiency  of 
optical-to-THz  power  conversion. 

B.  Work  Performed 

Work  performed  in  Phase  I  of  this  project  has  been  mostly  focused  on  theoretical  analysis  of  the 
proposed  optical  to  THz  energy  converter.  The  main  goal  has  been  to  suggest  a  configuration 
which  would  provide  conversion  efficiency  much  better  than  has  been  recently  reported.1 1 

We  have  analytically  and  numerically  solved  the  following  tasks: 

1.  Formulated  the  confinement  theory  for  nonlinear  transformation  and  proposed  multi¬ 
layer  waveguide  with  TEM  THz  wave,  which  maximizes  the  difference  frequency 
generation. 
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2.  Derived  analytically  dispersion  equations  for  THz  and  optical  modes. 

3.  Proved  that  phase-matching  condition,  which  is  necessary  for  efficient  energy  transfer 
from  optical  to  THz  wave,  is  satisfied  for  the  proposed  multi-layer  waveguide. 

4.  Proved  that  absorption  losses  for  TEM  THz  mode  are  in  acceptable  range. 

5.  Designed  a  grating  In-coupler  for  optical  beams  with  high  excitation  efficiency. 

6.  Derived  analytically  an  expression  for  energy  transformation  coefficient  and  found  its 
numerical  estimate  for  the  proposed  structure. 

C.  Results  Obtained 

Difference  frequency  generation  (DFG)  presents  a  clear  path  in  the  development  of  a  coherent 
THz  source.12  9]  In  this  technique  two  laser  beams  with  close  frequencies  are  mixed  in  a  second- 
order  nonlinear  crystal  to  produce  a  difference  THz  frequency.  A  typical  setup  widely  used  for 
wavelength  conversion  is  shown  in  Figure  l.[9]  It  has  a  pump  Nd:YAG  laser,  two  KTP  crystals 
placed  inside  a  resonator  (built  by  two  mirrors  Ml  and  M2)  for  two  close  optical  wavelength 
generation  and  a  bulk  BNA  crystal,  taken  as  an  example  of  a  nonlinear  crystal,  for  Difference- 
Frequency  Generation.  As  it  is  reported  in  [9],  these  authors  managed  to  achieve  13.3% 
efficiency  of  pump  to  idler  waves  conversion.  But  idler  to  THz  wave  conversion  has  still  been  a 
low  (-0.0001-0.001%)  efficiency  process.  [2~9] 

A  problem  with  low  conversion  efficiency  comes  from  the  weak  nonlinear  coupling  between 
optical  beams  and  THz  wave.  As  we  show  in  the  Figure  1  the  only  part  needed  to  be  improved  in 
the  current,  state  of  the  art,  DFG  THz  generation  method  is  the  way  of  excitation  of  a  nonlinear 
layer. 


where  a  change  is  proposed 


THz-wave 
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To  propose  an  efficient  optical-to-THz  converter  with  optimized  coupling  between  e/m  waves 
one  needs  to  formulate  an  analytical  theory  of  the  general  energy  transformation,  which  we 
successfully  delivered  in  this  project. 

Cl.  Analytical  theory  of  nonlinear  mixing  in  waveguides 

We  have  derived  electromagnetic  power  of  a  THz  wave  and  a  transformation  coefficient  in  an 
arbitrary  waveguide  configuration: 

Pthz  =  V  Paj1Paj2  (1) 


V  = 


87T2rie// _ 1_ 

n1n2nTHz  ce0AjHz  Aeff 


L2  sine2 


(2) 


where 


Aeff  ~ 


A* 


THz 


y(2)  =i  =i 

1^(2)  |  ■  Po)1Po)2  d x  dy  > 


(3) 


//  Ev  ■  dx  dy  =  SV[l,  (4) 

is  the  waveguide  mode  normalization;  L  is  the  length  of  energy  conversion  along  the  waveguide; 
P0)l,  P(I>2  are  the  electromagnetic  powers  of  the  optical  idler  modes  in  the  waveguide;  x('2)  is  the 
second-order  optical  susceptibility  tensor,  and  dejf  is  its  maximum  element  -  nonlinear  optical 
coefficient;  n1,n2,nTHz  are  the  effective  propagation  indexes  of  the  interacting  modes;  Akz  is 
the  wavevector  mismatch  between  the  three  modes;  Aeff  is  the  effective  area  which  quantifies 
the  magnitude  of  the  coupling  between  the  nonlinear  polarization,  excited  THz  mode  and  the 
optical  modes. 

Possible  absorption  losses  of  electromagnetic  waves  in  the  waveguide  were  omitted  in  Eq.2. 
Generally  they  are  required  to  be  considered.  In  this  report  we  lead  to  a  configuration  where 
absorption  losses  are  negligible. 

If  we  suppose  that  mode  configuration  of  the  waveguide  is  chosen  in  such  way  that  vector 
multiplication  produces  maximum  value,  what  will  be  achieved  by  correct  choice  of  the  material 
and  orientation-control  crystal  growth  (discussed  further  in  C.5),  then  the  problem  of  maximizing 
transformation  coefficient  splits  into  two  tasks: 

1.  Maximizing  overlap  integral  0  —  ff  ETHz  ECJiEa)2  dx  dy  with  the  normalization 
condition  for  interacting  modes  //  Ev  ■  E^  dx  dy  =  5V[l , 

2.  Providing  phase  matching  between  modes:  AkzL  — >  0, 

3.  Correcting  Overlap  equation  with  possible  absorption  losses  of  the  DFG  mode;  finding 
absorption  loss  for  proposed  configuration. 

Maximizing  the  overlap  between  modes: 

In  Figure  2  we  show  the  graphical  concept  of  finding  a  maximum  overlap  between  optical  and 
THz  modes.  It  is  clear  that  maximum  overlap  integral  achieved  in  case  of  fundamental  modes. 
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Electric  field  of  a  fundamental  mode  is  of  same  sign  in  the  whole  cross-section  of  a  waveguide. 
However,  low-order  optical  can  also  be  considered. 


Figure  2.  A  simple  confinement  model  for  overlap  between  optical  and  THz  modes  (under  the 
normalization  condition) 


For  a  planar  waveguiding  structure  we  obtain  the  following  equation  for  modes  overlap: 

n  V  2 

/  nx  \ 

( -  -  COS 


2opt 


0{CLTHz  >  aopt )  —  (  fa, 


"'J  aTHz  \aTHzJ  aopt 


■cos 


(— )]  dx 

\aopt/ _ 


maximum 


(5) 


When  deriving  Equation  5  we  assumed  that  both  idler  optical  beams  only  one  (same)  optical 
mode  in  a  nonlinear  optical  layer  of  the  waveguide,  which  is  true  because  their  frequencies  are 
very  close. 


Figure  3.  Overlap  at  aopt=1.5pm  (A,opt  ~1.5pm)  as  a  function  of  THz  wave  thickness 
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In  Figure  3  we  show  the  overlap  between  modes  at  1.5pm  optical  thickness  layer.  It  drastically 
decays  with  the  increase  of  THz  wave  thickness.  Ideally,  one  has  to  squeeze  THz  mode  into  a 
size  of  the  optical  wavelength  to  maximize  energy  conversion  between  them,  which  is  hard  to 
achieve. 

Another  interesting  and  most  surprising  result  we  get  from  mode  overlap  analysis  when  we  fix 
the  thickness  of  a  THz  layer.  In  Figure  4  we  show  such  case.  If  ,  suppose,  100  micron  thickness 
is  the  thinnest  which  we  are  to  realize  experimentally  for  terahertz  wave  then  there  is  no  reason 
to  make  the  optical  layer  very  thin  (down  to  1  micron),  overlap  value  does  not  increase  much. 


Figure  4.  Overlap  at  aTHz=  1 00p  m  as  a  function  of  the  optical  layer  thickness 


Following  these  results  we  can  state  that  the  original  idea  of  using  a  single  metallic  plate  covered 
with  multi-layer  dielectric  suggested  in  Phase  I  proposal  (Figure  5)  will  not  work  at  THz 
frequency  range  (0.1-10THz).  The  matter  is  in  the  confinement  of  the  surface  plasmon  polaritons 
mode  on  smooth  metals.  The  confinement  of  a  THz  surface  plasmon  in  perpendicular  to  metallic 
surface  direction  is  in  order  of  10cm,  which  is  very  week.  Overlap  integral  of  THz  wave  with 
optical  modes  would  be  very  small.  This  configuration  will  only  work  for  IR  (k=2-10micron) 
difference  frequency  generation. 


t  on  w'l'Mon  I  ninth  L 


Figure  5.  TE  and  TM  configuration  for  THz  generation  proposed  in  Phase  I  proposal.  They  will  only  work  for  b=2- 
lOmicron  difference  frequency  generation. 

Fortunately,  we  can  use  a  quasi-TEM  mode  which  exists  between  two  parallel  metallic  plates 
(Fig. 6)  at  any  spacing  between  them.[1H]  In  the  proposed  configuration  optical  beams  are 
confined  in  a  thin  film  of  nonlinear  dielectric.  There  are  two  concerns:  phase-matching  between 
THz  and  optical  modes,  and  the  absorption  losses  of  the  TEM  THz  mode.  As  it  was 
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experimentally  shown  in  several  papers[11'12'13'141  this  loss  is  very  low  ~  0.2dB/cm  for  lOOmicron 
spacing  between  metallic  plates.  Also  a  successful  confinement  of  THz  wave  with  frequency 
-ITHz  in  a  8micron  (!!)  thick  spacing  with  negligible  absorption  losses  has  been  reported.11’1 


Metallic  Layer 


Figure  6.  Proposed  TEM  mode  multi-layer  waveguide  for  optical  to  THz  converter. 


C2.  Dispersion  equation  for  TEM  THz  mode  and  satisfaction  of  the  phase¬ 
matching  condition 

It  is  important  to  derive  a  dispersion  equation  for  TEM  THz  mode  and  prove  that  phase¬ 
matching  between  THz  and  optical  modes  is  possible.  Also,  this  equation  will  allow  finding  the 
exact  absorption  losses  of  the  THz  signal. 

In  order  to  find  a  configuration  which  would  have  phase-matching  between  optical  and  proposed 
TEM  THz  mode  we  had  to  solve  a  dispersion  equation  for  5-layer  planar  symmetrical 
waveguide. 


TM  guided  modes  are  found  in  the  following  form: 


r  Ae~h3(x~D\ 

B cos(h2x)  +  T  sin(/t2x) , 
Hy(x)  —  \  A  cosh (/ijx)  +  E  sinhChjx) 
Z  cos(/t2x)  +  H  sin(/t2x) , 

geh3  (*+£>) 


x  >  D,n  —  n3 

d  <  x  <  D,n  —  n2  h\  —  /cf  —  koni 

—  d  <  x  <  d,n  —  %  ,  where  h\  —  — /cf  +  /cfnf 

—  D  <  x  <  —d,n  —  n2  hi  =  /cf  -  /cfnf 

x  <  —D,n  =  n3 


In  case  of  Quasi-TEM  mode  located  between  two  metallic  plates  we  need  to  put: 
n\  >  0,n|  >  0 ,Re(nl)  <  0 


We  find  the  dispersion  equation  for  TM  Even  modes: 


■jitanhCM)  =  -|{tan  [h2(D  -  d)  +  atan  (jj)  ]} 


(6) 


For  even  modes:  A  —  9 ,B  —  Z,T  —  —H,  E  —  0 

B  —  A  [cos (h2D)  +  ^  Jsin (7i2D)],  B  —  A  [sin (h2D)  -  ^^|cos(  h2D)j 


A  —  A 


cosh  (h-td) 


h3  nj 


cos(/i2(D  —  d))  +  - — 2  sin  h2(D  —  d) 


h2  n\ 
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1  /A‘ 


=  —( - - - .' 

n\  \cosh(/i xd), 


ho  rio 


-,2 


+  2 

ni 


D  —  d  ( 


cos(d2(D  —  d))  +  - — ^  sin(^2 (P  ~  d)) 

sin(2/t2(D  —  d)) 


h2  n\ 


hi  nf 


2  V 


I  1  +  ““4 


^2 


+ 


2d|  n| 


+ 


2  do  n: 


cos(2/i2(D  —  d))  — 


H  n\ 
4d|  nf 


4d2 

sin(2d2(D  —  d)) 


sinh(2/i1d)  d 
+  2 


+ 


2  Mi 


We  have  numerically  coded  Equation  8  for  our  configuration  A  1/S  i  CWGa A  s/S  i  Oi/ A 1 .  Parameters 
for  phase  matching  condition  between  optical  and  THz  modes  are  shown  in  the  Proposal. 


We  have  also  obtained  TEM  modes  by  FDTD  coding  (Fig. 7, 8). 


Hy 


Air  /  waveguide  boundary 


Figure  7.  Typical  picture  of  TEM  mode  at  A=300p  m.  (space  dimensions  are  in  microns,  amplitude  of  the 
fields  are  shown) 


Ex 


Figure  8.  Typical  picture  of  TEM  mode  at  A.=10pm,  (space  dimensions  are  in  microns,  real  value  of  the 
fields  are  shown) 
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C3.  Found  phase-matched  configuration 

For  difference  frequency  TFIz  generation  PMC  can  be  written  as  following: 

<Winz(6)i)  -  w2nz(n>2)  =  ruTHzn.z (n>THz)  , 

where  co12  are  the  frequencies  of  the  two  optical  idler  beams,  coTHz  —  a)1  —  0)2  is  the  difference 
frequency,  nz(m)  is  the  effective  refraction  index  of  the  propagating  mode.  In  general,  the 
structure  of  all  three  modes  can  differ  very  much.  However,  in  our  case  —  nz(m2)  - 

both  idler  beams  excites  the  same  mode  in  the  waveguide.  So,  the  PMC  looks  simpler: 

nz(o)  i)  =  nz(mTHz)  (7) 

As  it  turns  out  satisfaction  of  PMC  is  feasible  for  both  proposed  configuration.  Our  optical 
beams  are  confined  in  the  Layer  2  (Fig. 6).  After  getting  a  THz  TM  mode  dispersion  equation  for 
5  layer  system:  metal/dielectric/nonlinear-dielectric/dielectric/metal,  optical  TE  mode  dispersion 
equation  for  3  layer  system  dielectric/nonlinear-dielectric/dielectric  and  numerically  solving 
them  we  have  obtained  the  following  configuration  that  exhibit  a  phase-matching  between  THz 
and  optical  modes: 

Al(thickness  10pm)/Si02(4.8pm)/GaAs(0.385pm)/Si02(4.8pm)/Al(10pm). 

At  optical  excitation  wavelengths  A1  —  1.5pm,  X2  =  1.5  —  1.515pm  there  is  a  full  match 
between  effective  refractive  indexes  in  this  optical  and  0.5-3THz  range  (nen=2. 15).  In  Figure  9 
we  have  plotted  these  modes:  TEM  (l=300pm),  and  TEi  (1=1. 5pm). 


TE\1 


Ref  Index,  Ey  Optical 


Position  [um] 


Figure  9.  TEM  THz  mode  (l=300pm),  and  Optical  (1=1. 5pm)  mode  in  the  found  phase-matched 
configuration 
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C4.  Absorption  losses  of  the  TEM  THz  wave 

k 

To  find  the  absorption  losses  we  seek  —  =  nTHz  =  Re(nTHz)  +  i/m(nTHz). 

For  example,  let  us  find  the  absorption  loss  for  the  found  phase-matched  configuration: 
Al(thickness  10pm)/Si02(4.8pm)/GaAs(0.385pm)/Si02(4.8pm)/Al(10pm).  At  optical  excitation 
wavelengths  At  —  1.5  p.m,  X2  —  1.5  —  1.51 5 p.m  there  is  a  full  match  between  effective 
refractive  indexes  in  this  optical  and  0.5-3THz  range  ( Re(nTHz )  =  neJJtical  —  2.146  ). 

We  find  the  “TEM  THz”  root  of  the  Dispersion  Equation  6  graphically.  We  have  prepared  a 
numerical  code  for  this,  and  in  the  Figure  10  we  show  how  the  imaginary  part  of  nTHz  is  found. 
At  A/thz  =214pm  we  get  nTHz  =  2.146+i0.0073. 


Figure  10.  Finding  “TEM  THz”  root  graphically  k=214pm.  x=Re(nTH/). 

Now  we  can  calculate  the  absorption  losses: 

^hz  =  -El  exp(— t/c0ttTHzz)  ~ exp(aTHz) , aTHz  =  /c0/m(nTHz)  (8) 

For  the  considered  case  (A^214pm)  we  get  the  conversion  length  Lc  =  aTHz_1  =  4.7mm,  which 
is  a  remarkable  huge  length  for  optical  effects!  Therefore,  the  proposed  coherent  energy 
converter  can  have  up  to  1cm  waveguide  length,  which  is  achievable  by  CVD  and  ALD  growing 
methods. 
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C5.  Final  design  of  microstructured  multi-layer  waveguide  and  grating  In¬ 
couplers 

Based  on  the  results  described  above  we  have  formulated  two  possible  design  of  the  proposed 
optical-to-THz  coherent  energy  converter.  As  it  turns  out  two  configurations  are  possible:  TE 
optical  and  TM  optical  excitations. 

TE  optical  excitation: 

As  it  shown  in  the  Figure  1 1  the  Electric  field  Ey  of  the  TE  optical  modes  excites  polarization 

->  NL  — 

Px  of  a  43m  point  symmetry  group  nonlinear  semiconductor  (GaAs,  GaP),  which  is,  in  turn, 
parallel  to  the  electric  field  of  TEM  difference  frequency  wave.  Grating  coupler  will  be  etched 
into  the  SiOa  layer  (Layer  1).  Upper  metal  layer  is  not  “covering”  grating  grooves,  so  the  laser 
light  can  efficiently  penetrate  Si02  layer  and  excite  the  grating. 


Figure  11.  TE  optical  -  TEM  THz  e/m  energy  converter  based  on  43m  point  symmetry  group  nonlinear 
material  (GaAs,  GaP,  etc.) 

Special  requirements  have  to  be  met  for  GaAs  layer. [15]  It  has  to  be  (100)  grown  layer  with  a 
grating  grooves  parallel  to  (110)  plane  of  the  crystal.  In  Figure  12  we  show  crystal  orientation  for 
the  case  of  TE  optical  excitation.  Fortunately,  there  is  a  developed  technique  to  grow  (100)  GaAs 
layer.[lf>l  In  Ref.l  ITO  coated  (200nm  thickness)  quartz  substrates  were  used  for  synthesis  of  thin 
film  GaAs. 

Grating  In-coupler 


Figure  12.  The  orientation  of  43m  crystal  for  converting  TE  optical  input  into  TEM  THz  wave.  Top 
surface  is  (100) 
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TM  optical  excitation; 

In  this  case  (Fig.  13)  the  dominant  electric  field  in  the  TM  wave,  Ex  (Ex  »  Ez  )  excites  the 

-»  NL 

nonlinear  polarization  Px  of  3m  symmetry  group  material  (LiNb03,  LiTa03,  etc).  LiNb03 
needs  to  be  special  c-oriented  sputter  deposited  on  SiOi.  Fortunately,  this  technology  is  already 
well  developed.117’18’19’20’211 


Figure  13.  TM  optical  -  TEM  THz  converter  based  on  3m  symmetry  group  nonlinear  material  (LiNhOs, 
LiTaCF,  etc) 

Grating  In-couper; 

Before  we  estimate  the  conversion  efficiency  for  the  proposed  optical-to-THz  converter  we 
should  provide  a  design  for  a  grating  coupler.  Grating  coupler  serves  as  an  efficient  exciter  of  the 
optical  modes  in  the  thin  film  nonlinear  layer.  The  PI  of  this  project  has  already  delivered 
numerically  and  experimentally  such  an  efficient  coupler.  The  measured  efficiency  was  80%.  By 
other  words,  80%  of  incident  optical  power  was  “inserted”  into  a  thin  (325nm  thickness)  optical 
waveguide  (A,=671nm).  1221 

Here  we  need  to  prove  numerically  that  we  can  have  this  high  efficiency  couplers  for  any 
wavelength,  particularly  at  A„pt=l  .5pm,  proposed  as  a  pumping  wavelength  for  GaAs  nonlinear 
layer. 

In  Figure  14  we  show  a  very  promising  design  of  such  a  coupler  which  works  at  Aop(=l  .5pm,  and 
have  coupling  efficiency  o  60%  at  the  normal  (!)  incidence.  The  period  of  the  grating  is  700nm, 
thickness  of  the  GaAs  waveguiding  layer  is  385nm,  and  groves  depth  is  120nm.  Groves  width  is 
350nm.  As  of  now  a  6  pm  y-wide  beam  simulated.  However,  as  PI  already  proved 
experimentally/221  the  efficiency  stays  high  even  for  very  y-wide  beam  (up  to  several  hundred 
microns). 
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Poynting  Vector^ 


11.25 


7.50 


3.75 


15.00 


Figure  14.  Grating  coupler  simulated  by  FDTD  numerical  method  at  X=  1.5pm.  x,z  are  in  microns.  As 
clearly  seen,  the  optical  mode  TE!  is  being  excited  apparently  only  in  +z-direction  at  normal  (!)  incidence. 


C6.  Energy  transformation  coefficient  in  general  case  with  absorption  losses 

We  have  derived  a  general  equation  for  energy  transformation  efficiency.  Here  we  show  a  brief 
derivation  of  this  equation. 

In  undepleted  case,  when  optical  powers  P  1,2  do  not  change  much  along  the  propagation 
(waveguide)  direction  we  have: 
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rf^THz 

dz 


—  ~  “thz^thz  +  kA1A*2e 


-iAk7z 


(9) 


where  ATHz  is  the  amplitude  of  the  THz  wave,  A1 2,  are  the  amplitudes  of  two  optical  modes,  k  is 


proportional  to  overlap  between  modes,  k 


ff* 


THz 


y(  2)  =i  =, 

_± _ .  u  u 

\X(2)\  ■  C<"1C"2 


dx  dy 


,  Akz  is  the 


wave  vector  mismatch  between  optical  and  THz  modes,  is  aTHz  the  absorption  loss  of  the  THz 
mode  mainly  occurring  in  the  metallic  layers  (Fig. 6). 

Solution  of  the  Differential  Equation  1  we  seek  in  the  following  form: 

^thzO)  =  /(z)e~a™zZ, 

which  lead  to: 

^  =  /cM;e(a™z-iAfcz)z  . 


dz 


After  integrating  (9)  we  obtain  the  equation  for  energy  transformation  coefficient: 


T]  = 


8n2djff 


(l+e~ 


£THzL-2e' 


cos(A/czL)) 


pip2  n^riTHz  C£0; If  Hz  Aeff 


(aTHz2+Akz  ) 


(10) 


where  L  is  the  mode  interaction  length;  Aejj- 


\x(2)^:  E(o1koj2 


dx  dy 


|  JT  E'VHz 

overlap  surface  between  modes,  which  obeys  the  mode  normalization  condition: 


is  the 


//  Ev(x,y)  ■  (x,y)  dx  dy  =  5V(1, 

n1;  n2,  nTHz  are  the  effective  propagation  indexes  of  the  interacting  modes;  /(2>  is  the  second- 
order  optical  susceptibility  tensor,  and  deJf  is  its  maximum  element  -  nonlinear  optical 
coefficient. 


Equation  10  is  very  general  equation  that  describes  coherent  energy  transfer  between  two  optical 
idler  modes  and  the  THz  mode  in  any  waveguide  geometry.  The  only  assumption  of  undepleted 
regime  for  optical  modes  can  be  corrected  in  our  further  research. 


This  equation  takes  into  account  both  factors  which  limit  the  efficiency  of  the  energy 
transformation:  absorption  losses  and  the  phase  mismatch  between  the  modes.  As  we  have 
shown  in  C.3  of  this  report,  we  have  found  a  configuration  which  provides  the  phase-matching 
condition:  A kzL  ->  0.  In  this  case  the  only  limiting  factor  is  the  absorption  loss  of  the  THz  mode. 
Equation  10  will  transform  into  the  following: 


v  =  h^=  ^2d2eff  i _ i_ 

/  pip2  nin2nTHz  c^o^thz  Aeff  (<*thz2) 


(11) 


Here  we  assumed  that  the  length  of  the  waveguide  will  be  slightly  more  than  the  conversion 

1 

length,  which  in  the  phase-matched  case  equal  to  Lc  — - . 

®THz 
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Applying  Equation  11  to  the  proposed  TEM  THz  waveguide: 

It  would  be  of  great  help  for  this  project  to  have  an  analytical  expression  for  the  energy 
transformation  coefficient  in  terms  of  parameters  of  the  proposed  waveguide  design  (Fig. 6).  A 
unique  property  of  the  TEM  mode,  which  is  its  step-function  dependence  on  the  normal 
coordinate  (x-axis,  Fig.  6,9),  helped  us  derive  such  an  equation. 


By  introducing  the  characteristic  optical  beams  space  confinement  in  y-direction  Ly,  we  have 


i 

Aeff 


JT Ethz  (x,y’)E(0l(x,y)E(iy2(x,y)  dx  dy\  =  L2y  |/J  ETHz  (x)FWi(x)£^2(x)  dx 


(12) 


Recalling  that  the  optical  frequencies  are  very  close:  a)1  -  <n2  =  &)THz  «  aq  2,  we  have  their  mode 
structures  the  same:  EWi(x)  =  E0>2  (x). 


Optical  modes  are  confined  in  a  thin  nonlinear  layer,  where  TEM  THz  electric  field  is  constant. 
So,  we  get: 

J^}-l2y(^ L(*  =  *Gaj)  dx. 

After  applying  the  normalization  condition  TEM  THz  mode  we  obtain: 

M~f  =  (^™z  =  *GaA^)  ’ 

and  using  D  field  boundary  condition: 


D 


THZ 

X,  GaAs 


=  D 


THZ 
X,  Si02 


we  finally  obtain: 


1 

Aeff 


Energy  conversion  efficiency  will  be  found  as 


1 


GaAs 


T]  = 


PiP- 


1^2 


8  n2d2eff 

n1n2nTHz  C£o^-thz 


1 

(aTHz2)  ' 


(13) 


Now  we  can  easily  predict  the  conversion  efficiency  coefficient  for  our  multi-layer  waveguide 
once  we  know  the  exact  absorption  losses. 

Estimate  of  energy  conversion  efficiency: 

We  can  finally  calculate  the  expected  energy  conversion  efficiency.  We  will  assume  that  grating 
coupler  provides  maximum  -100%  excitation  efficiency.  Eq.  13  transforms  into  following: 
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2d 


\ 


V  = 


eff 


PtP2  n^riTHz  c£q 


Lyh  GaAs 


1  + 


(K^Y(hSi  qA 

vn™2;  UGaJ 


/ 


(Im(nTHz)2) 


(14) 


It  is  interesting  that  Eq.14  does  not  depend  explicitly  on  optical  idler  beam  wavelengths,  neither 
on  generated  THz  wavelength.  So,  we  obtain: 

For  LiNb03  at  Ai ,2=1 ,5|jm 


2- ( 175- 10  12) 


n  = 


3- 108  8.8510  12  2-2-2-0.00^-10G10  6  0.385 10  6 


1  + 


5.2  )  '  (l0- 10  6  -0.38510  6) 


2.11 


,0. 


38510  6) 


=  3.249x  10 


-6 


For  GaAs  at  Ai ,2=1 ,5|jm 


2- (50- 10  12) 


n  = 


3- 108  8.85 10  12  2-2-2-0.0052-10010  6  0.385 10 


,-6 


1  + 


3.6 

2.11 


,10  10  6  -  0.385 10  6) 


,0. 


38510 


■-6) 


=  1.15x  10 


Obtained  conversion  efficiency  numbers  are  not  optimized.  We  think  that  by  squeezing  optical 
beam  in  y-direction  down  to  30  micron  and  additional  squeeze  of  TEM  mode  down  to  1-2 
microns  in  x-direction  can  improve  efficiency  up  to  10  4  value. 

Also,  there  are  organic  DAST  and  BNA  crystals  which  have  deff  ~103pm/V.|23  24’25]  If  they  could 
be  employed  in  our  micron  thick  waveguide  geometry  -  they  would  increase  q  up  to  10  2  value. 
The  best  value  of  conversion  efficiency  recently  reported  in  Ref.  [1]  is  1.3- 10"  W"  . 

D.  Estimate  of  Technical  Feasibility 

In  this  Phase  I  SBIR  project  we  have  performed  a  vigorous  theoretical  analysis  of  the  proposed 
multi-layer  waveguide.  To  our  best  knowledge  no  one  yet  in  the  scientific  world  has  done  any 
similar  analysis.  This  research  effort  outlines  and  finds  the  best  design  configuration  for 
Difference  Frequency  Generation  (DFG). 

The  objective  of  this  particular  project  has  been  to  formulate  a  design  for  THz  DFG.  But  the 
approach  which  we  have  developed  within  this  project  is  not  limited  to  DFG  and  neither  to  THz 
frequency  range  as  a  goal.  This  approach  will  be  the  best  for  deep  UV,  IR  (FWIR,  MWIR, 
SWIR)  radiation  generation.  It  will  lead  to  some  unique  waveguiding  designs  for  these  e/m  wave 
ranges.  The  designs  will  be  different  from  the  one  for  THz  presented  here,  but  they  will  be  also 
the  most  efficient. 
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In  this  report  we  have  not  specified  any  sources  of  the  idler  optical  beams.  The  choice  of  them  is 
arbitrary.  It  is  clear  that  the  best  performance  will  be  with  coherent  laser  beams,  since  they  can 
be  squeezed  easily. 

Another  important  aspect  is  the  coherence  of  the  THz  signal  generated  in  the  proposed  converter. 
Time  coherence  of  the  THz  signal  will  be  transferred  (translated)  from  idler  optical  signal 
without  any  change.  There  is  no  electronic  transition  between  bands  involved  in  this  nonlinear 
conversion  technique  as  contrary  to  “photoconductive  antenna”  method[26]  where  electron 
scattering  limits  coherence  of  the  DFG  signal.  In  our  scheme  the  narrower  the  optical  lasing 
wavelength,  the  narrower  the  generated  THz  line. 

We  have  already  mentioned  that  the  success  of  this  project  now  lies  in  advanced  area  of  the 
microfabrication.  We  have  specified  the  references  where  thin  LiNb03  or  GaAs  layers  were 
successfully  grown  and  tested  on  their  values  of  nonlinear  coefficient. 

This  project  is  concerned  with  development  of  tunable  time-coherent  THz  source  of  0.1-10  THz 
frequency  range.  In  the  further  fabrication  stage  of  this  research  two  types  of  coherent  THz 
source  will  be  developed:  a  quasi-cw  tunable  narrowband  source  with  frequency  line  width 
100MHz  (line  quality  factor  ~104,  repetition  rate  ~10kHz),  and  a  cw  broadband  tunable  source 
with  frequency  line  width  ~0.3THz.  These  two  lasers  will  be  based  on  different  pump  lasers: 

OPO  with  a  solid  state  single  mode  laser  and  with  two  DFB  lasers,  which  are  frequency  tuned  by 
temperature.  Tunability  range  of  will  be  0.1/0.3THz  -  10  THz.  Our  unique  multi-layer  waveguide 
will  provide  us  the  maximum  possible  efficiency  of  the  nonlinear  electromagnetic  energy 
conversion  from  optical  to  THz  frequency  domains.  We  expect  to  reach  an  average  output  power 
of  lmW  required  by  the  solicitation. 

The  successful  completion  of  this  project  in  the  future  will  benefit  many  new  commercial 
markets.  Having  coherent,  tunable  and  handy  THz  source  will  bring  scientists  an  opportunity  to 
fully  explore  THz  band.  Particular  applications  include:  spectroscopic  miniature  sensors  of  THz 
spectra  of  chem/bio-moleculas;  next  generation  THz  electronics  components  with  aim  to  future 
THz  networks  and  computing  chips;  vision  through  clothing  at  the  airport  for  security  purposes; 
non-destructive  imaging  through  teeth,  since  THz  easily  penetrates  them;  biological 
hyperspectral  imaging  of  human  body. 
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